arXivil509.01687vl [cond-mat.mes-hall] 5 Sep 2015 


Single Molecule Transistor based Nanopore for the detection of 

Nicotine 

S. J. Ray* 

Institute of Materials Seienee, Teehnical University of Darmstadt, 
Alarich-Weiss-Str. 2, 64.287 Darmstadt, Germany 

Abstract 

A nanopore based detection methodology was proposed and investigated for the detection of 
Nicotine. This technique uses a Single Molecular Transistor (SMT) working as a nanopore op¬ 
erational in the Coulomb Blockade regime. When the Nicotine molecule is pulled through the 
nanopore area surrounded by the Source(S), Drain(D) and Gate electrodes, the charge stability 
diagram can detect the presence of the molecule and is unique for a specific molecular structure. 
Due to the weak coupling between the different electrodes which is set by the nanopore size, the 
molecular energy states stay almost unaffected by the electrostatic environment that can be re¬ 
alised from the charge stability diagram. Identification of different orientation and position of the 
Nicotine molecule within the nanopore area can be made from specific regions of overlap between 
different charge states on the stability diagram that could be used as an electronic fingerprint for 
detection. This method could be advantageous and useful to detect the presence of Nicotine in 
smoke which is usually performed using chemical chromatography techniques. 
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I. INTRODUCTION 


Nicotine is a colourless organic compound that is often considered toxic [1] when present at 
a significant concentration. Environmental tobacco smoke (ETS) resulting from the burning 
of a cigarette contains Nicotine [2] and the exposure to this is considered to be an impor¬ 
tant reason behind several chronicle diseases like lung cancer [3-5], asthma [6-8] and heart 
diseases [9-11]. When exposed to air, the presence of Nicotine stays significantly strong for 
around 2 hours and gets metabolised to Cotine at a later stage. The detection of Nico¬ 
tine is thus a sensitive technique as ETS is the only possible source of Nicotine in air [2]. 
Usually the detection of Nicotine is performed using the techniques of radioimmunoassay 
[12, 13], liquid chromatography [14], gas chromatography [15] etc. Due to the advance¬ 
ment in the spectrographic and chromatographic techniques, the mostly used methods use a 
combination of gas chromatography, mass spectrometry and high performance liquid chro¬ 
matography with flame ionisation or nitrogen specific detectors [2]. These techniques are 
sensitive to the amount of Nicotine and can detect a smaller presence of Nicotine in air. How¬ 
ever, they usually require a series of chemical preparations and analysis and sophisticated 
instrumentations. 

Here a simple but effective method was proposed for the detection of Nicotine using 
a nanopore based sensing technique. Conventionally a nanometer sized hole is created 
artificially via exposing ion beam [16] or electron beam [17] on an electrically insulating 
membrane to create a solid state nanopore. Due to the smaller size of the nanopore, when a 
molecule passes through a nanopore connected between two electrodes in an ionic solution, 
a change of the current can be detected. This detection principle is commonly used for the 
sequencing of DNA [18]. The materials used for constructing a nanopore are varied like 
Silicon [17], Silicon Nitride [16], Protein [19] and Graphene [20] used in the recent times. 
However, the amplitude of the tunnelling current generated in the presence of a molecule 
within a nanopore depends on the coupling between them, which is usually small, thus 
limiting the detection efficiency. 

By way of exploiting the weak-coupling property, a nanopore based approach has been 
proposed for the detection of DNA [24]. Here in this work, the nanopore is essentially a 
single molecular transistor with the source/drain and gate electrodes placed on three sides 
of the nanopore. The Nicotine molecule is the active component of the device which is 
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equivalent to the ‘island’ in a single electron transistor. The operation of a SMT can be 
described using the Orthodox theory of Coulomb Blockade [21] which has been explained in 
details here [23]. The usefulness of a SMT has been demonstrated for charge detection [25] 
and tracking the the impurity location in a single molecular device [26] . 

In this case, the Nicotine molecule was placed at the centre of the nanopore and the 
nanopore was made in a way with one of the 3-dimensions is significantly smaller than the 
other two. Thus when the Nicotine molecule is pulled through the SMT, the presence of the 
molecule can be found from the charge stability diagram. Due to the sequential nature of 
transport operation, the SMT is sensitive to the orientation of the molecule and this adds 
additional sensitivity to the operation of the device. In this work, it was demonstrated that 
the charge stability diagram is unique for different orientation of the molecule and mildly 
sensitive to the position of the molecule within the SMT-nanopore which could be used as 
an electronic fingerprint for the detection of Nicotine. 


II. COMPUTATIONAL RECIPE 

The SMT based nanopore used in the present work has a rectangular structure which 
was illustrated in Fig. 1(b). The Source (S) and the Drain (D) electrodes were placed on 
top of a dielectric layer at its two ends. The dielectric layer has a thickness of 4 A with 
dielectric constant (e) = IOsq- A metallic gate electrode of thickness 1 A was placed at the 
back of the dielectric layer covering the entire area. All the electrodes were considered to be 
made of Gold with work function (IF) = 5.28 eV and can be replaced by other metals too. 
The S/D electrodes have a thickness of 5 A along the a:-axis and the horizontal separation 
(along the .s-axis) between them is 15 A. The height of the S/D electrodes along the r/-axis 
is 12 A and the nanopore has a pore cross section of 12 Ax 15 A area. Along the a:-axis, 
the entire structure has a thickness of 5 A which is comparable to the Silicon [16, 17] and 
Graphene[20] based nanopores reported earlier. 

Nicotine is chemically known as G10H14N2 with the pyridine and pyrrolidone rings at¬ 
tached to each others. This structure used for the present work was obtained using Open 
Babel package [22] as illustrated in Fig. 1(a). Due to the non-planar structure of the molecule 
with respect to the a:y;2-coordinates, the electrostatic environment is different for different 
orientations of the molecule within the nanopore area. Owing to the weak-coupling regime 
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FIG. 1: (a) Schematic of the Nicotine molecule used in the current work. The atoms in ‘grey’, 
‘white’ and ‘blue’ colours represent Carbon, Hydrogen and Nitrogen atoms respectively as shown 
in the labels, (b) A sample SMT device structure (^90) with the molecule placed at the centre 
and gate connected to the back of the dielectric layer, (c) Total energy of the molecule as function 
of the gate voltage for different charge states for the ^90 configuration and (d) Charge stability 
diagram for the same in ^90 configuration, the colorbar on the right represents different charge 
states. 

of operation, the molecular energy levels stay almost unaffected by the electrostatic envi¬ 
ronment. This is ensured by way of making a larger nanopore area which also allows the 
possibility of easier fabrication of such structure. Such a device is operational when the 
molecule is present within the nanopore area which can be realised from the charge stability 
diagram that is different from the situation when the molecule is absent within the SMT. 
In an experimental situation, the molecule could arrive within the nanopore in a variety of 
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FIG. 2: Charge stability diagram of the molecule within the nanopore-SMT for various configu¬ 
rations as follows (a) xO, (b) x90, (c) xl80, (d) x270, (e) 2/180, (f) 2/270, (g) z90, (h) 2:180 and 
(i) 2:270. Schematic of the nanopore with the molecule was illustrated in (a) for xO, (f) for 2/270 
and (i) for 2:270 configurations. The orientation of the molecule in a x 2 / 2 :-coordinate frame were 
shown for the other configurations in (b) x90, (c) xl80, (d) x270, (e) 2/180, (g) 2:90 and (h) 2:180 
configurations. 
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FIG. 3: (a) Charging energy obtained from the reference diamond (marked by ‘d’ in Fig. Id and 
Fig. 2) for different configurations of the molecule within the SMT, (b) Fitted parameters a and 
/3 (from Eqn.l) for these various configurations for the q=3 charge state. 


orientations of its atoms with respect to the electrodes and the corresponding response of 
the SMT will carry the signature of it. This was investigated systematically by way of rotat¬ 
ing the molecule along the 3 cartesian axes and finding the influence of different molecular 
configurations on the charge stability diagram of the SMT. 

Different configurations were named according to different rotational angles of the 
molecule with respect to the cartesian axes. Starting from an initial configuration (called 
as xO as illustrated in Fig. 2(a)) [29], these configurations were named following the angle 
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FIG. 4: (a) Vg^ref for different configurations of the molecule as estimated from the point ‘A’ in 
Fig. 2, (b) Angular dependence of Vg^ref for a continuous rotation of the molecule around the x-axis 
keeping the other coordinates fixed. 


of rotation from its initial state, some of which were listed in Table. I. The electrostatic 
behaviour of the SMT has been investigated by estimating the energy of the molecule using 
first-principle based calculations within a DFT framework. This computational approach 
was introduced by Kaajesberg et. al [28] and formulated within a DFT framework by Stok- 
bro [23]. The charging energy estimated by this method has been found to be in excellent 
agreement with the experimental results in various systems [23]. For the present work, 
Ab-initio calculations were performed by using Atomistic Toolkit Package [27] developed 
by the Quantum Wise division. The self consistent calculations were performed under the 


7 












































Configuration 

Description 

xO 

Default/initial configuration. 

x90 

Rotated by 90° along the x-axis from xO 

configuration keeping (^, z) coordinates fixed. 

xl80 

Rotated by 90° along the x-axis from x90 

configuration keeping (^, z) coordinates fixed. 

2/90 

Rotated by 90° along the y-ax.is from xO 

configuration keeping (x, z) coordinates fixed. 

z270 

Rotated by 270° along the 2 ;-axis from xO 

configuration keeping {x^y) coordinates fixed. 


TABLE I: Description of different device configurations 


Local density Approximation (LDA) within a non-spin polarised DFT framework where the 
wave functions were expanded within a double-!^ polarised (DZP) basis set. All the metallic 
electrodes were considered to be equipotential surfaces and Neumann Boundary Conditions 
were used for solving the Poisson’s equation. 

III. RESULTS AND DISCUSSION 

The total energy of the molecule in the SMT for the y90 configuration was plotted as 
function of the gate voltage in Fig. 1(c). For different charge states, the total energy contains 
a fixed part contributed by the metal electrodes {=eV) which for the present case is 5.28 eV. 
Under the influence of a positive gate voltage, negative charge states have a lower energy 
and the reverse happens for the positive charge states at negative gate voltages. For a 
specific charge state q, the total energy varies almost linearly with the gate voltage and the 
dependence can be further expressed analytically by, 

Eiq) = Eo + aqVg + /3{eVgf ( 1 ) 

The first term Eq is the zero-term in energy that is independent of the gate voltage, e 
is the electronic charge and the proportionality constants a and {3 could be estimated by 
fitting Eqn. 1 to the computationally estimated energies. For the present case, the average 
value of a has been found to be 0.3 and P ~ 0.0013 eV“^ from the fits as illustrated in 










Fig. 1(c). The value of a quantifies the strength of direct electrostatic coupling between 
the gate electrode and the molecule, while the second term is a measure of the overall 
polarisation of the different charge states under the electrostatic environment of the SMT. 
For a planar structure of the molecule lying parallel to the dielectric layer, the value of 
(3 is usually smaller than the present case as the molecular structure is not entirely two 
dimensional which induces additional polarisation for different charge states. However, the 
value of a ^ indicates the gate dominated linear regime of operation of the current device. 

The conduction behaviour of the SMT nanopore can be understood in details from the 
charge stability diagram as illustrated in Fig. 1(d) for the y90 configuration. Due to the 
sequential nature of transport in the CB regime, a finite conduction in such a device only 
occours for certain values of the applied Source-Drain Bias (Vd) that comes within the 
accessible molecular energy levels. In order to observe a finite non-zero conduction, the 
charge population of the island has to be changed by integer numbers which can be done by 
changing Vd (keeping Vg fixed) or changing Vg for a specific value of Vd- For a fixed value 
of Vd, an electron can be added or removed from the island by changing Vg. In the other 
way, at a specific value of Vg, the charge state of the island can be changed to an excited 
state via changing Vd- A series of these line scans for a certain range of the other parameter 
results in diamond shaped regions on the Vd — Vg plane that separates the non-conducting 
and conducting areas of the 2D map. The large central diamond marked by ‘d’ in Fig. 1(d) 
represents the neutral state of the molecule and the neighbouring diamonds represent the 
‘cationic’ and ‘anionic’ state of the molecule to its ‘left’ and ‘right’ respectively. The charging 
energy of the molecule in the ground state for the neutral case can be estimated from the 
height of the central diamond (marked by ‘d’) which is 8.88 eV for y90 configuration. 

However, the charge stability diagram does not stay the same when the position and 
orientation of the molecule changes within the SMT which was illustrated in details in 
Fig. 2. Starting from the xO configuration, the charge stability diagram changes for different 
configurations of the molecule when it is rotated along different Cartesian axes. System¬ 
atic evolution of the different charge states can be found from the charge stability diagram 
through the evolution of different coloured regions to follow the changes in different configu¬ 
rations. By looking at the charge stability diagrams for different configurations with naked 
eyes, several features could be noticed. The major changes between different configurations 
come from the relative orientation of the molecule (and the atoms) with respect to the elec- 
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trodes which is reflected in the charging energy (Ec) behaviour of the molecule in various 
configurations. This can be compared for different configurations by considering a specific 
diamond on the charge stability diagram which for the present case has been marked by ‘d’ 
in Fig. 2. Systematic comparison of the charging energy estimated from the height of the 
diamond ‘d’ has been illustrated in Fig. 3(a). 

It is to be noted that the energy levels on the charge stability diagram represents the 
charging energy of the molecule, which is the difference between two energy states [= E{n + 
1) — E(n)] and is independent of the electrical polarisation contribution towards the energy. 
Thus the variability of Ec does not contain any non-linear energy term from Eqn. 1 and can 
be related to the effect of direct gate-molecule coupling for various configurations. From the 
dependence of the total energy with the gate voltage for the charge state q =3, the values 
of the fitted parameters a and /3 for various configurations were plotted in Fig. 3(b). For 
different configurations, the overall variation of (5 is much smaller compared to the change 
in a and the variation of Ec can be compared to the variation in a. 

For a given configuration of the molecule within the SMT, a increases when a larger 
volume of the molecule (or the atoms) stays in closer proximity of the dielectric layer covering 
the gate electrode. The difference in the values of Ec (in Fig. 3(a)) between the xO and the 
a:90 configurations can be related to a similar change in the value of a in Fig. 3(b). By 
rotating the molecule further by 90° from the x90 configuration, the shape and size of ‘d’ 
changes with an increase in Ec- Due to the symmetry of the molecular structure with respect 
to the electrodes, Ec for the x9 and xl80 configurations are of similar values as for the a:90 
and x270 cases which can be found in Fig. 3(a). This can be compared to the value in 
a in Fig. 3(b) which also goes through similar variations between the xO, a:90, xl80, x270 
configurations. By making similar comparisons, Ec for the y90 and y270 were estimated to 
be the same as for the cases of xO and ?/180. However, the situation is not the same for 
various configurations while rotating the molecule around 2 :-axis. The Ec for ^90 is different 
from z270 case, as in one case, the Hydrogen atoms attached to the hexagonal side faces the 
electrode while for the other case, the pentagonal arm comes close to the dielectric layer. 
The influences of gate coupling a is higher for z90 as a larger numbers of atoms come in 
close proximity to the gate electrode than the ^270 case. In Fig. 3(b), the value of a for 
xO, xl80, r/180, z90 are of similar orders which is the same for the case of Ec for these cases 
in Fig. 3(a). Between the {y90,y270) and (a;0,r/180) configurations, a larger variation in a 
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can be observed compared to the variations for other configurations. 

In this nanopore based detection approach, the Nicotine molecule is pulled through the 
SMT-nanopore and the orientation of the molecule can be different from the present situ¬ 
ation. This possibility was considered for the present case and a detection procedure was 
proposed that could consider these various situations. As illustrated in Fig. 2, the overall 
charge stability diagram changes for different configurations and specific regions of them 
could be considered to understand the effect of various orientations. In order to quantify 
this evolution in a systematic manner, a reference point was considered on the charge stabil¬ 
ity diagram which was marked by point ‘A’ in Fig. 2. This reference point was chosen to be 
the left vertex of the extreme left diamond where charge states q = 0 and q = 1 intersects 
each other. The a:-coordinate of ‘A’ (i.e. the value of Vg) has been plotted in Fig. 4(a) to 
understand the evolution of the reference point for different configurations. The differences 
in the values of I4,ref for different orientations of the molecule around a:-axis is relatively 
smaller than the differences in I4,ref for similar orientations around the y and 2: axes. The 
value of 14,ref for the xO and xl80 configurations are of similar order as it is the same for a:90 
and a:270 configurations. The contrast in I^.ref increases significantly between the y90 and 
yl80 cases. For different configurations around the 2-axis, the value of I^.ref can be found 
to be different for the 290, 2I8O and 2270 cases. 

Fig. 4(a) considers only a few discrete configurations and to study further the dependence 
of the charge stability diagram for various rotational angles the molecule, systematic study 
was performed for configurations with different rotation angles around the a:-axis keeping 
the other two coordinates fixed. The values of I4,ref estimated for these cases were plotted as 
function of different rotation angles in Fig. 4(b). The value of Ig,ref does not overlap between 
two neighbouring configurations. Although a discrete set of data points were considered for 
this case, but the trend indicates that the value of I4,ref is periodic for 180'^ rotation. As 
the value of I4,ref is unique for different configurations, thus this could be considered as 
a method of identification of different configurations of the molecule within the nanopore. 
Since it is easy to obtain the charge stability diagram experimentally by measuring the 
Source-Drain current (Ids) for different values of the Vd and Vg, this could be used as an 
electronic fingerprint for the detection of Nicotine. As the SMT-nanopore works in the weak 
coupling limit, so the features in the charge stability diagrams are mainly determined by the 
energy state of the molecule. Therefore the colour pattern and the charge stability diagram 
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carry intrinsic signatures of the presence of Nicotine molecule in a SMT-nanopore uniquely. 

As Nicotine is the strongest component present in the smoke, hence a direct exposure 
of the smoke could be detected using this technique efficiently. However, for the case of 
experimental or practical realisation of such a method, several factors need to be considered. 
The present case only considers the situation when only Nicotine molecule is present within 
the SMT-nanopore area which is not necessarily the situation as the smoke also contains 
other elements of air. However, the charge stability diagram is unique for each molecule and 
hence the presence of Nicotine could be easily identified with background subtraction by this 
method. In the present case, the device performance was estimated at 300K and the effect 
of thermal noise will be present at this temperature of operation. However, since the voltage 
range of operation is significantly higher (V^ ~ 1 Volt), this contribution could be ignored. 
There is another possibility of fluctuation arising from the movement of the molecule itself 
within the nanopore. This is going to add a much smaller influence on the charge stability 
diagram as the weak-coupling operation of this device will have a minor influence due to 
this. 


IV. CONCLUSION 

In this work, a novel approach based on a SMT based nanopore was proposed and inves¬ 
tigated for the detection of Nicotine. The first principal based calculation was performed on 
such a system to obtain the charge stability diagram in the CB regime. Due to weak cou¬ 
pling regime of operation of this SMT-nanopore, the charge stability diagram carries unique 
signatures of the molecule which could be used as an electronic fingerprint for the detection 
of Nicotine. For various orientations of the molecule, the charge stability diagram is dis¬ 
tinctive for different configurations that adds additional sensitivity for identifying specific 
configurations of Nicotine within the nanopore area. This device becomes operational when 
the molecule is present within the nanopore area and the corresponding charge stability 
diagram is easy to obtain experimentally. Further investigation was performed by focusing 
on certain specific regions of the charge stability diagrams, where different charge states 
meet and the value of I4,ref for these positions were found to be unique for each different 
configurations of the molecule within the nanopore. This methodology could work as a very 
powerful tool for the detection of Nicotine and Nicotine related compounds present in the 
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smoke. 
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